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ABSTRACT 

We present an analysis of Spitzer-IRS spectroscopic maps of the LI 157 pro- 
tostellar outflow in the H2 pure- rotational lines from S(0) to S(7). The aim of 
this work is to derive the physical conditions pertaining to the warm molecular 
gas and study their variations within the flow. The mid-IR H 2 emission follows 
the morphology of the precessing flow, with peaks correlated with individual CO 
clumps and H2 2.12/im ro- vibrational emission. More diffuse emission delineat- 
ing the CO cavities is detected only in the low-laying transitions, with Ji ower < 
2. The H 2 line images have been used to construct 2D maps of N(H 2 ), H 2 
ortho-to-para ratio and temperature spectral index /3, in the assumption of a gas 
temperature stratification where the H 2 column density varies as T _/3 . Variations 
of these parameters are observed along the flow. In particular, the ortho-to-para 
ratio ranges from ~0.6 to 2.8, highlighting the presence of regions subject to 
recent shocks where the ortho-to-para ratio has not had time yet to reach the 
equilibrium value. Near-IR spectroscopic data on ro- vibrational H 2 emission have 
been combined with the mid-IR data and used to derive additional shock param- 
eters in the brightest blue- and red-shifted emission knots. A high abundance of 
atomic hydrogen (H/H 2 ~ 0.1-0.3) is implied by the observed H 2 column densities, 
assuming n(H 2 ) values as derived by independent SiO observations. The pres- 
ence of a high fraction of atomic hydrogen, indicates that a partially-dissociative 
shock component should be considered for the H 2 excitation in these localized 
regions. However, planar shock models, either of C- or J-type, are not able to 
consistently reproduce all the physical parameters derived from our analysis of 
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the H 2 emission. Globally, H 2 emission contributes to about 50% of the total 
shock radiated energy in the L1157 outflow. We find that the momentum flux 
through the shocks derived from the radiated luminosity is comparable to the 
thrust of the associated molecular outflow, supporting the scenario where this 
latter is driven by the shock working surface. 

Subject headings: STARS: FORMATION, ISM: JETS AND OUTFLOWS, ISM: 
INDIVIDUAL OBJECTS: L1157, INFRARED: ISM: LINES AND BANDS 

1. Introduction 

The interaction of protostellar outflows with the ambient molecular cloud occurs through 
radiative shocks that compress and heat the gas, which in turn cools down through line 
emission at different wavelengths. In the dense medium where the still very embedded 
protostars (the so called class sources) are located, shocks are primarily non-dissociative, 
and hence the cooling is mainly through emission from abundant molecules. Molecular 
hydrogen is by far the most abundant species in these environments, and although H2 emits 
only through quadrupole transitions with low radiative rates, it represents the main gas 
coolant in flows from young protostars. H 2 shocked emission in outflows has been widely 
studied in the past mainly through its ro-vibrational emission in the near-IR (e.g. Eisloffel 
et al. 2000; Giannini et al. 2004; Caratti o Garatti et al. 2006) that traces the dense gas at 
T~2000-4000 K. Most of the thermal energy associated with the shocks is however radiated 
away through the emission of H 2 rotational transitions of the ground state vibrational level 
at A < 28/im (e.g. Kaufman & Neufeld 1996). Mid-IR H 2 lines are easily excited at low 
densities and temperatures between 300 and 1500 K: therefore they are very good tracers of 
the molecular shocks associated with the acceleration of ambient gas by matter ejection from 
the protostar. Given the low excitation temperature, they can also probe regions where H 2 
has not yet reached the ortho-to-para equilibrium, thus giving information on the thermal 
history of the shocked gas (Neufeld et al. 1998; Wilgenbus et al. 2000). In addition, given 
the different excitation temperature and critical densities of the v=0-0 and v> 1 H 2 lines, 
the combination of mid-IR with near-IR observations is a very powerful tool to constrain the 
global physical structure and the shock conditions giving rise to the observed emission. 

The study of the 0-0 rotational emission in outflows started in some detail with the 
Infrared Space Observatory. Thanks to the observations performed with the SWS and ISO- 
CAM instruments, the shock conditions, the ortho-para ratio and the global H 2 cooling have 
been derived in a handful of flows (e.g. Neufeld et al. 1998; Nisini et al. 2000; Molinari et 
al. 2000; Lefloch et al. 2003). More recently, the Infrared Spectrometer (Houck et al. 2004) 
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on board Spitzer, with its enhanced spatial resolution and sensitivity with respect to the 
ISO spectrometers, has been used to obtain detailed images of the H 2 rotational emission, 
from S(0) to S(7), of several outflows, from which maps of important physical parameters, 
such as temperature, column density and o/p ratio, have been constructed (Neufeld et al. 
2006; Maret et al. 2009; Dionatos et al. 2010). In this framework, Neufeld et al. (2009, 
hereafter N09) have recently presented IRS spectroscopic maps observations of five young 
protostellar outflows at wavelengths between 5.2 and 37/xm, and discussed their averaged 
physical properties and overall energetics. In all the flows, the H 2 S(0)-S(7) emission has 
been detected and contributes to more than 95% of the total line luminosity in the 5.2-37yum 
range, while atomic emission, in the form of Fell and SI fine structure lines, accounts for 
only the remaining ~5%. 

In the present paper, we will analyse the H 2 line maps obtained by N09 towards the 
LI 157 outflow, with the aim of deriving the main physical conditions pertaining to the 
molecular gas and their variations within the flow. This in turn will give information on the 
thermal history of the flow and on how energy is progressively transferred from the primary 
ejection event to the slow moving ambient gas. 

For this first detailed analysis, L1157 has been chosen among the sources observed by 
N09 given its uniqueness as a very active and well studied flow at different wavelengths. 
More than 20 different chemical species have been indeed detected in the shocked spots 
of this object (Bachiller & Perez Gutierrez 1997; Benedettini et al. 2007), some of them 
for the first time in outflows (e.g. HNCO, Rodriguez- Fernandez et al. 2010, and complex 
organic molecules, Arce et al. 2008, Codella et al. 2009) , testifying for a rich shock induced 
chemistry. Warm H 2 shocked emission in LI 157 is also evidenced through near-IR maps 
(e.g. Davis & Eisloffel, 1995) and Spitzer-IRAC images (Looney et al. 2007). The LI 157 
outflow has been also recently investigated with the Herschel Space Observatory, showing to 
be very strong also at far-IR wavelengths (Codella et al. 2010, Lefloch et al. 2010, Nisini et 
al. 2010). 

The L1157 outflow extends about 0.7 pc in length. Its distance is uncertain and has been 
estimated between 250 and 440 pc. Here we will adopt D=440 pc for an easier comparison 
with other works. The outflow is driven by a highly embedded, low mass class source 
(L1157-mm or IRAS20386+6751) having L boi ~ 8.3 L & (Froebrich 2005). It is a very nice 
example of an outflow driven by a precessing and pulsed jet, possessing an S-shaped structure 
and different cavities, whose morphology has been reproduced assuming that the outflow is 
inclined by ~80° to the line of sight and the axis of the underlying jet precesses on a cone of 
6° opening angle (Gueth et al. 1996). The episodic mass ejection events are evidenced by the 
presence, along the flow, of individual clumps that are symmetrically displaced with respect 
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to the central source. It is therefore a very interesting target for a study of the physical 
conditions pertaining to these active regions through an H 2 excitation analysis. 

The paper is organized as follow: the observations and the main results are summarized 
in §2. In §3 we describe the analysis performed on the H 2 images to derive maps of tempera- 
ture, column density and ortho-to-para ratio. A more detailed NLTE analysis on individual 
emission peaks is also presented here, where the Spitzer data are combined with near-IR 
data to further constrain the excitation conditions. The implications of these results for the 
shock conditions along the LI 157 flow are discussed in §4, together with an analysis of the 
global energy budget in the flow. A brief summary follows in §5. 

2. Observations and results 

Observations of the LI 157 outflow were obtained in November 2007 with the IRS in- 
strument, during Cycle 4 of the Spitzer mission . The full IRS spectral range (5.2-36.5/im) 
was observed with the Long-High (LH), Short-High (SH) (R ~ 600) and Short-Low (SL) (R 
between 64 and 128) modules. The LI 157 outflow region was covered through 5 individual 
IRS maps of ~ l'xl' of size each, arranged along the outflow axis. Each map was obtained 
by stepping the IRS slit by half of its width in the direction perpendicular to the slit length. 
For the SH and LH modules the slit was stepped also parallel to its axis by 4/5 (SL) and 
1/5 (LH) of its length. Details on the data reduction that generated the individual line 
maps from the IRS scans are given in N09. The final maps have been resampled to a grid 
of 2" spacing allowing a pixel by pixel comparison of maps obtained with the different IRS 
modules. Maps of the brightest detected lines as well as the full spectrum in a representative 
position are shown in Fig. 7 and Fig. 12 of N09. As regards to H 2 , all the pure rotational 
lines of the first vibrational levels, from S(0) to S(7), are detected at various intensity along 
the flow. Here we report, in Tab. [H the H 2 brightness measured in a 20" FWHM Gaussian 
aperture towards different positions. 

Fig 1 shows the L1157 maps of the S(l) and S(2) lines while Fig. 2 displays the S(5) line 
with superimposed contours of the CO 2-1 emission from Bachiller et al. 2001. In the same 
figure, a map of the 2.12 /zm 1-0 S(l) line is also presented. The morphology of the 0-0 S(5) 
and 1-0 S(l) is very similar, with peaks of mid-IR emission located at the near-IR knots 
from A to D, as identified by Davis & Eisloffel (1995). When compared with the CO map, 
the mid-IR H 2 emission appears to follow the curved chain of clumps (labelled as B0-B1-B2 
and R0-R1-R, for the blue-shifted and red-shifted lobes, respectively) that also correspond 
to peaks of SiO emission, as resolved in interferometric observation by Gueth et al. (1998) 
and Zhang et al. (2000). The LI 157 outflow morphology has been suggested to delineate 
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a precessing flow (Gueth et al. 1996), where the H 2 and SiO peak emission knots follow 
the location of the actual working surface of the precessing jet and are thus associated with 
the youngest ejection episodes. Diffuse H 2 emission is also detected in the S(l)-S(2) maps, 
that delineates the wall of a cavity that connects the central source with both the RO and 
BO clumps. Such a cavity has been recognized in the CO 1-0 interferometric maps and it is 
likely created by the propagation of large bow-shocks. The S(l)-S(2) maps of Fig. 1 show 
extended emission of H 2 also in the SE direction (i.e. where the B2 clump is located) and 
in the eastern edge of the northern lobe, that also follow quite closely the CO morphology: 
these regions at lower excitation might trace additional cavities created by an older ejection 
episode of the precessing jet. 



3. H 2 Analysis 

3.1. LTE 2D analysis of the rotational lines: maps of averaged parameters 

We have used the H 2 line maps to obtain the 2D distribution of basic H 2 physical 
parameters, through the analysis of the rotational diagrams in each individual pixel. As 
described in N09, who analysed the global H 2 excitation conditions in LI 157, the distribution 
of upper level column densities of the S(0)-S(7) lines as a function of their excitation energy, 
does not follow a straight line, indicating that a temperature stratification in the observed 
medium exists. The exact form of this temperature stratification depends on the type of 
shock the H 2 lines are tracing, as they probe the post-shock regions where the gas cools 
from ~ 1000 K to ~ 200 K. 

The simplest way to parametrize the post-shock temperature stratification is to assume 
a power- law distribution: this approach was applied by Neufeld & Yuan (2008), which also 
show that this type of distribution is expected in gas ahead of unresolved bow-shocks. On 
this basis, and following also N09, we fit the observations assuming a slab of gas where the 
H 2 column density in each layer at a given T varies as dN oc T~@dT. 

This law is integrated, to find the total column density, between a minimum (T min ) 
and a maximum (T max ) temperatures . For our calculations we have kept T max fixed at 
4000 K, since gas at temperatures larger than this value is not expected to contribute to 
the emission of the observed pure rotational lines. T m i n was instead assumed to be equal, 
in each position, to the minimum temperature probed by the observed lines. This T m i n is 
taken as the excitation temperature giving rise to the observed ratio of the S(0) and S(l) 
column densities, assuming a Boltzman distribution. The T min value ranges between ~150 
and 400 K. 
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We found that the approach of a variable T min produces always fits with a better x 2 than 
assuming a fixed low value in all positions. We also assume the gas is in LTE conditions. 
Critical densities of rotational lines from S(0) to S(7) range between 4.9 cm -3 (S(0)) and 
4.4xl0 5 cm" 3 (S(7)) at T=1000 K assuming only H 2 collisions (Le Bourlot et al. 1999): 
critical densities decrease if collisions with H and He are not negligible. Deviations from 
LTE can be therefore expected only for the high-J S(6) and S(7) transitions: the S/N of 
these transitions in the individual pixels is however not high enough to disentangle, in the 
rotational diagrams, NLTE effects from the effects caused by the variations of the other 
considered parameters. In particular, as also discussed in N09, there is a certain degree of 
degeneracy in the density and the /3 parameter of the temperature power law in a NLTE 
treatment that we are not able to remove in the analysis of the individual pixels. This issue 
will be further discussed in §3.21 An additional parameter of our fit is the ortho-to-para ratio 
(OPR) value. It is indeed recognized that the 0-0 H 2 lines are often far from being in ortho- 
to-para equilibrium, an effect that in a rotational diagram is evidenced by a characteristic 
zigzag behavior in which column densities of lines with even- J lie systematically above those 
of odd- J lines. In order not to introduce too many parameters, we assume a single OPR 
value as a free parameter for the fit. In reality, the OPR value is temperature dependent 
(e.g. N09 and §3.2[) . and therefore the high-J lines might present an OPR value closer to 
equilibrium then the low- J transitions. Our fit gives therefore only a value averaged over 
the temperature range probed by the lines considered (i.e. ~200-1500 K). 

In summary, we have varied only three parameters, namely the total H 2 column density 
N(H 2 ), the OPR, and the temperature power law index (3, in order to obtain the best model 
fit through a x 2 minimization procedure and assuming a 20% flux uncertainty for all the 
lines. The fit was performed only in those pixels where at least four lines with an S/N 
larger than 3 have been detected. Before performing the fit, the line column densities were 
corrected for extinction assuming A v = 2 (Caratti o Garatti et al. 2006) and adopting the 
Rieke & Lebofsky (1985) extinction law. At the considered wavelengths, variations of A v of 
the order of l-2mag do not affect any of the derived results. 

Figure [3] shows the derived map of the H 2 column density, while in Fig. H] and [5] maps 
of the OPR and are displayed together with temperature maps relative to the "cold" and 
"warm" components (T coW and T warm ), i.e. the temperature derived from linearly fitting 
the S(0)-S(l)-S(2) and S(5)-S(6)-S(7) lines, once corrected for the derived OPR value. In 
Fig. 6 we also show the individual excitation diagrams for selected positions along the flow, 
obtained from intensities measured in a 20" FWHM Gaussian aperture centered towards 
emission peaks (Tab. 1). Values of the fitted parameters in these positions are reported in 
Tab. [2j In addition to T co u and T warm , we give in this table also the values of the average 
temperature in each knot, derived through a linear fit through all the H 2 lines (T me( j). 
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The maps show significant variations in the inferred parameters along the outflow. The 
H 2 column density ranges between 5xl0 19 and 3xl0 20 cm -3 . The region at the highest 
column density is located towards the Bl molecular bullet (see Fig. [2]o. This is consistent 
with the higher column density of CO found in Bl with respect to other positions in the 
blue lobe (Bachiller & Perez Gutierrez 1997), and might suggest that this is a zone where 
the outflowing gas is compressed due to the impact with a region of higher density (Nisini 
et al. 2007). Towards the NW, red-shifted outflow, the column density has a more uniform 
distribution, with a plateau at ~ 10 20 cm -2 that follows the H 2 intensity distribution. The 
N(H2) decreases at the apex of the red-shifted outflow, with a value slightly below 10 20 cm -2 
at the position of the D near-IR knot. 

Tcoia ranges between ~ 250 and 550 K. The highest values are found at the tip of the 
northern outflow lobe, while local maxima corresponds to the positions of line intensity 
peaks. T warm ranges between ~ 1000 and 1500 K. In this case the highest values are in 
the southern lobe, at the position of the A NIR knot. As a general trend, the T warm value 
decreases going from the southern to the northern peaks of emission, with the minimum 
value at the position of the D NIR knot. 

/3 values range between ~4-4.5 in the blue-shifted lobe while it is larger in the red-shifted 
lobe, with maximum values of ~ 5.5 at the tip of the flow. Due to the degeneracy between (3 
and density discussed in the previous section, these values can be considered as upper limits 
because of our assumption of LTE conditions. Neufeld & Yuan (2008) have discussed the j3 
index expectations in bow-shock excitation. A /3 index of ~ 3.8 is expected in paraboloid 
bow shocks having a velocity at the bow apex high enough to dissociate H2, in which case the 
temperature distribution extends to the maximum allowed temperature. Slower shocks that 
are not able to attain the maximum temperature, produce steeper temperature distributions, 
i.e. with values of /3 greater than 3.8. This is consistent with our findings: low values of (3 
(of the order of 4) are found in the blue-shifted lobe: here, evidence of H 2 dissociation is 
given by the detection of atomic lines (i.e. [Sill], [Fell] and [SI]) in the IRS spectrum. In 
the red-shifted lobe, where values of /3 larger than 4 are derived in the LTE assumption, 
no atomic emission is detected and the T warm values are lower than those measured in the 
blue-shifted lobe, indicating a maximum temperature lower than in the blue flow . 

The OPR varies significantly along the outflow, spanning from ~ 0.6 to 2.8. Hence 
it is always below the equilibrium value of 3. Although a one-to-one correlation between 
temperature and OPR cannot be discerned, some trends can be inferred from inspection of 



1 Fig. [5] shows that the molecular clumps Bl, R0 and R coincide in position with the NIR H2 knots A, C 
and D. In the paper, both nomenclature will be used specifying if we refer to the NIR or mm condensations 



-8 - 



Fig. [5j In general the OPR minima are observed in plateau regions between two consecutive 
intensity peaks, where also the cold temperature has its lowest values. At variance with 
this trend, the emission filaments delineating the outflow cavity within ± 20" from the mm 
source, where the cold temperature reaches a minimum value of ~ 250 K, show rather high 
values of the OPR, ~2.4-2.8. This might suggests that this region has experienced an older 
shock event that has raised the OPR, though not to the equilibrium value, and where the 
gas had time to cool at a temperature close to the pre-shock gas temperature. On the other 
hand, at the apex of both the blue- and red-shifted lobes, where the cold temperatures are 
relatively high (i.e. 500-550 K), the OPR is rather low, 1.5 — 2.0. Evidence of regions of 
low OPR and high temperatures at the outflow tips has been already given in other flows 
(Neufeld et al. 2006; Maret et al. 2009). It has been suggested that these represent zones 
subject to recent shocks where the OPR has not had time yet to reach the equilibrium value. 



3.2. NLTE analysis: constraints on H and H 2 particle density 

Additional constraints on the physical conditions responsible for the H2 excitation, 
are provided by combining the emission of the mid-IR H2 pure-rotational lines from the 
ground vibrational level with the emission from near-IR ro-vibrational lines. It can be seen 
from Fig. [2] that the 2.12/xm emission follows quite closely the emission of the 0-0 lines 
at higher J. In addition to the 2.12/im data presented in Fig. [2J we have also considered 
the NIR long-slit spectra obtained on the A and C NIR knots by Caratti o Garatti (2006). 
These knots, at the spatial resolution of the NIR observations (i.e. ~0.8"), are separated 
in several different sub-structures that have been individually investigated with the long-slit 
spectroscopic observations. For our analysis we have considered the data obtained on the 
brightest of the sub-structures, that coincide, in position, with peaks of the 1-0 S(l) line. 

In order to inter-calibrate in flux the Spitzer data with these NIR long-slit data, obtained 
with a slit-width of 0.5 arcsec, we have proceeded as follows: we first convolved the 2.12/im 
image at the resolution of the Spitzer images and than performed photometry on the A 
and C peaks positions with a 20" diameter Gaussian aperture, i.e. with the same aperture 
adopted for the brightness given in Tab. [TJ We have then scaled the fluxes of the individual 
lines given in Caratti o Garatti (2006) in order to match the 2.12/im flux gathered in the slit 
with that measured by the image photometry. In doing this, we assumed that the average 
excitation conditions within the 20" aperture are not very different from those of the A-C 
peaks. This assumption is observationally supported by the fact that the ratios of different 
H2 NIR lines do not change significantly (i.e. less than 20%) in the A-C knot substructures 
separately investigated in Caratti o Garatti (2006). We have considered only those lines 
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detected with S/N larger than 5; in practice this means considering lines from the first 
four and three vibrational levels for knots A and C, respectively. The excitation diagrams 
obtained by combining the Spitzer and NIR data for these two knots are displayed in FigJHJ 

In order to model together the 0-0 lines and the near-IR ro-vibrational lines, we have 
implemented two modifications to the approach adopted previously. First of all, the NIR 
lines probe gas at temperatures higher than the pure rotational lines, of few thousands of 
K, at which values it is expected that the OPR has already reached equilibrium. Thus, the 
ortho-para conversion time as a function of temperature needs to be included in the fitting 
procedure, since lines excited at different temperatures have different OPRs. We have here 
adopted the approach of N09 and used an analytical expression for the OPR as a function of 
the temperature, considering a gas that had an initial value of the ortho-to-para ratio OPRo 
and has been heated to a temperature T for a time r. Assuming that the para-to-ortho 
conversion occurs through reactive collisions with atomic hydrogen, we have: 



OPR(t) _ OPRp n (H)kr , OPRlTE /j _ n (H)kr\ Q\ 

l + OPR(r) 1 + OPR 1 + OPR L te [ ' 

In this expression, n(H) is the number density of atomic hydrogen and OPR^te is the 
ortho-to-para ratio equlibrium value. The parameter k is given by the sum of the rates 
coefficients for para-to-ortho conversion (k po ), estimated as 8xlO _11 exp(-3900/T) cm 3 s _1 , 
and for ortho-to-para conversion, k op ~ k po /3 (Schofield et al. 1967). Thus the dependence of 
the OPR on the temperature is implicitly given by the dependence on T of the k coefficient. 
The inclusion of a function of OPR on T, introduces one additional parameter to our fit: 
while we have previously considered only the average OPR of the 0-0 lines, we will fit here 
the initial OPRo value and the coefficient K = n(H)r. 

The second important change that we have introduced with respect to the previous 
fitting procedure, is to include a NLTE treatment of the H 2 level column densities. In fact, 
the critical densities of the NIR lines are much higher than those of the pure rotational lines 
(see Le Bourlot et al. 1999). For example, the n cr i t of the 1-0 S(l) 2.12/zm line is 10 7 cm -3 
assuming only collisions with H 2 and T=2000 K. Therefore the previously adopted LTE 
approximation might not be valid when combining lines from different vibrational levels. 
This is illustrated in FigJTJ where we plot the results obtained by varying the H2 density 
between 10 3 and 10 7 cm~ 3 while keeping the other model parameters fixed . The observed 
column densities in the A position are displayed for comparison. For the NLTE statistical 
equilibrium computation we have adopted the H 2 collisional rate coefficients given by Le 
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Bourlot et al. (1999) □ . This figure demonstrates the sensitivity of the relative ratios 
between 0-0 and 1-0 transitions to density variations. For example, in this particular case, 
the ratio N(H 2 ) _ 0S(7) /N(H 2 ) 1 _ 0S{1) is 64.6 at n(H 2 ) = 10 4 cm" 3 and 1.9 at n(H 2 ) > 10 7 
cm -3 . The figure also shows that the observational points display only a small misalignment 
in column densities between the 0-0 lines and the 1-0 lines, already indicating that the 
ro-vibrational lines are close to LTE conditions at high density. 

In FigJHl we show the final best-fit models for the combined mid- and near-IR column 
densities in the A and C positions. As anticipated, the derived n(H 2 ) densities are large, of 
the order of 10 7 and 6xl0 6 cm -3 , for the A and C positions, respectively. The two positions 
indeed show very similar excitation conditions: only the column density is a factor of 3 
smaller in knot C. Hence, we conclude that the lack of detection of rotational lines with 
v> 3 in knot C, in contrast to knot A (Caratti o Garatti et al. 2006), is due purely to a 
smaller number of emitting molecules along the line of sight and not to different excitation 
conditions. 

The derived H 2 densities are much higher than previous estimates based on other tracers. 
Nisini et al. (2007) derive a density of 4xl0 5 cm -3 at the position of knot A from multi-line 
SiO observations, thus more than an order of magnitude smaller than those inferred from our 
analysis. The high-J CO lines observed along the blue-shifted lobe of L1157 by Hirano et al. 
(2001), indicate a density even smaller, of the order of 4xl0 4 cm -3 . SiO is synthesized and 
excited in a post-shock cooling zone where the maximum compression is reached, therefore 
it should trace post-shock regions at densities higher than H 2 (e.g. Gusdorf et al. 2008). 
One possibility at the origin of the discrepancy is our assumption of collisions with only H 2 
molecules, and thus of a negligible abundance of H. This can be considered roughly true in 
the case of non- dissociative C-shocks, where H atoms are produced primarely in the chemical 
reactions that form H 2 from O and H 2 , with an abundance n(H)/n(H 2 +H) ~ 10~ 3 (e.g. 
Kaufmann & Neufeld 1996). However, if the shock is partially dissociative, the abundance 
of H can increase considerably and collisions with atomic hydrogen cannot be neglected, in 
view of its large efficiency in the H 2 excitation. This situation cannot be excluded at least 
for the knot A, where atomic emission from [Fell] and [SI] has been detected in our Spitzer 
observations. Since we cannot introduce the n(H) as an additional independent parameter of 
our fit, we have fixed n(H 2 ) at the value derived from SiO observations (4xl0 5 cm -3 , Nisini 
et al. 2007) and varied the H/H 2 abundance ratio. The best fit is in this case obtained with 
a ratio H/H 2 =0.3: this indicates that our observational data are consistent with previous H 2 
density determinations only if a large fraction of the gas is in atomic form. 



2 The rate coefficients for collisions with ortho- and para-H2, HI and He, computed in Le Bourlot et al. 
(1999) are available at the web-site: |http: / / ccp7. dur . ac . uk/cooling_by _h2/| 
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Turning back to the inferred OPR variations with temperature, our fit implies that the 
OPR in the cold gas component at T=300 K is significantly below the equilibrium value, 
while the value of OPR=3 is reached in the hot gas at T=2000 K traced by the NIR lines. 
The parameter K = n(H)r is constrained to be ~10 6 and 10 7 yrcm -3 for knots A and C, 
respectively. We can also estimate the time needed for the gas to reach this distribution of 
OPR, from the limits on the atomic hydrogen abundance previously discussed. Our data 
implies a high value of the n(H) density: a minimum value of n(H) ~0.6-lxl0 4 cm~ 3 , (for 
knots C and A, respectively) is given if we assume H/H 2 ~ 10 -3 (and thus the n(H 2 ) ~ 
10 7 cm -3 , given by our fit), while a maximum value of ~10 5 cm~ 3 , is derived from the fit 
where n(H 2 ) is kept equal to 4xl0 5 cm -3 . The high abundance of atomic H ensures that 
conversion of para- to ortho-H2 proceeds very rapidly: the fitted values of the K parameter 
indeed imply that the observed range of OPR as function of temperature have been attained 
in a timescale between 100 and 1000 yrs for both the knots. 

Finally, given the column density and particle density discussed above, we can estimate 
the H 2 cooling length (L ~ N(H 2 )/n(H 2 )). If we consider the case of n(H 2 ) ~ 10 7 cm~ 3 , and 
negligible n(H), we have L ~ 10 13 cm while a length of ~ 10 15 cm is inferred in the case of 
n(H 2 )~ 4xl0 5 cm -3 . All the parameters derived from the above analysis are summarised in 
Tab. [3] and they will be discussed in the next section in the framework of different shock 
models. 



4. Discussion 

4.1. Shock conditions giving rise to the H 2 emission 

The copious H 2 emission at low excitation observed along the LI 157 outflow indicates 
that the interaction of the flow with the ambient medium occurs prevalently through non- 
dissociative shocks. Both the Spitzer IRS maps of N09, and the NIR narrow band images of 
Caratti o Garatti et al. 2006, show that significant gas dissociation in LI 157 occurs only at 
the A peak, where both mid-IR lines from [Fell], [SII] and [SI] and weak [Fell] at 1.64/zm 
have been detected. Weaker [SI] 25/zm and [Fell] 26/im emission have been also detected on 
the C spot, but overall the atomic transitions give a negligible contribution to the total gas 
cooling, as pointed out in N09. These considerations suggest that most of the shocks along 
the outflow occur at speeds below ~ 40km s" 1 , as this is the velocity limit above which H 2 is 
expected to be dissociated. The knot A is the only one showing a clear bow-shock structure. 
Here the velocity at the bow apex is probably high, causing H 2 dissociation and atomic line 
excitation, consistent with the fitted temperature power law (3 of ~ 4, as discussed in §3.1, 
while the bulk of the H 2 emission comes from shocks at lower velocities originating in the 
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bow wings. 

Constraints on the shock velocity that gives rise to the molecular emission in L1157 have 
been already given in previous works. The sub-mm SiO emission and abundances, measured 
in different outflow spots, suggest shock velocities of the order of 20-30kms _1 (Nisini et al. 
2007). The comparison of SiO and H2 emission against detailed shock models performed by 
Gusdorf et al. (2008) confirm a similar range of velocities in the NIR-A knot, although the 
authors could not find a unique shock model that well represents both the emissions. 

Cabrit et al. (1999) found that the column density of the mid-IR H2 emission lines, from 
S(2) to S(8), observed by ISO-ISOCAM was consistent either with C-shocks having velocities 
of ~25kms~ 1 or with J-shocks at lower velocity, of the order of ~ 10km s _1 . Gusdorf et al. 
(2008), however, conclude that stationary shock models, either of C- or J-type, are not able 
to reproduce the observed rotational diagram on the NIR-A position, constructed combining 
ISOCAM data and NIR vibrational lines emission. A better fit was obtained by these authors 
considering non-stationary shock models, which have developed a magnetic precursor but 
which retain a J-type discontinuity (the so-called CJ shocks, Flower et al. 2003). Similar 
conclusions, but on a different outflow, have been reached by Giannini et al. 2006 who 
studied the H2 mid- and near-IR emission in HH54: in general, stady-state C- and J-type 
shocks fail to reproduce simultaneously the column densities of both the ro- vibrational and 
the v=0, pure-rotational H2 levels. 

A different way to look at the issue of the prevailing shock conditions in the observed 
regions, is to compare the set of physical parameters that we have inferred from our analysis 
to those expected from different shocks. With this aim, we summarize in Tab. [3] the physical 
properties derived on the A and C H2 knots. In addition to the parameters derived from 
the NLTE analysis reported in Section [21 namely H 2 post-shock density, H/H 2 fraction, 
initial OPR, cooling length and time, the table reports also the average values of OPR 
and rotational temperature, as they are measured from a simple linear fit of the rotational 
diagrams presented in Fig. [6j 

As mentioned in §3.2, the high fraction of atomic hydrogen inferred by our analysis rule 
out excitation in a pure C-shock. In fact, dissociation in C-shocks is always too low to have 
a H/H.2 ratio higher than 5xl0~ 3 , irrespective from the shock velocity and magnetic field 
strength (Kaufman & Neufeld 1996; Wilgenbus et al. 2000). C-shocks are not consistent with 
the derived parameters even if we consider the model fit with the high H 2 post-shock density 
of the order of 10 7 cm _3 and negligible atomic hydrogen: in this case we derive an emission 
length of 10 13 cm, which is much lower than the cooling length expected in C-shocks, which, 
although decreasing with the pre-shock density, is never less than 10 15 cm (Neufeld et al. 
2006) . 
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Stationary J-shock models better reproduce some of our derived parameters. For exam- 
ple, in J-type shocks the fraction of hydrogen in the post-shocked gas can reach the values 
of 0.1-0.3 we have inferred, provided that the shock velocity is larger than ~ 20 kms -1 . 
In general, a reasonable agreement with the inferred post-shock density and H/H 2 ratio is 
achieved with models having t> s =20-25 kms -1 and pre-shock densities of 10 3 cm -3 (Wilgenbus 
et al. 2000). Such models predict a shock flow time of the order of 100 yr or less, which 
is also in agreement with the value estimated in our analysis at least in knot A. In such 
models, however, the cooling length is an order of magnitude smaller than the inferred value 
of ~10 15 cm. In addition, the gas temperature remains high for most of the post-shocked 
region: the average rotational temperature of the v=0 vibrational level is predicted to be, 
according to the Wilgenbus et al. (2000) grid of models, always about 1600 K or larger, as 
compared with the value of about 800-900 K inferred from observations. The consequence of 
the above inconsistencies is that J-type shocks tend to underestimate the column densities 
of the lowest H 2 rotational levels in LI 157, an effect already pointed out by Gusdorf et al. 
(2008). 

As mentioned before, Gusdorf et al. (2008) conclude that the H 2 pure rotational emission 
in LI 157 is better fitted with a non-stationary C+J shock model with either v s between 20 
and 25 kms -1 and pre-shock densities tin = 10 4 cm -3 , or with v s ~ 15 kms -1 and higher 
pre-shock densities of nu = 10 5 cm -3 . Such models, however, still underestimate the column 
densities of the near-IR transitions: the post-shocked H 2 gas density remains lower than the 
NIR transitions critical density and the atomic hydrogen produced from H 2 dissociation is 
not high enough to populate the vibrational levels to equilibrium conditions. 

The difficulty of finding a suitable single model that reproduce the derived physical con- 
ditions is likely related to possible geometrical effects and to the fact that multiple shocks 
with different velocities might be present along the line of sight. It would be indeed inter- 
esting to explore whether bow-shock models might be able to predict the averaged physical 
characteristics along the line of sight that we infer from our analysis. 

4.2. Flow energetics 

H 2 emission represents one of the main contributor to the energy radiated away in shocks 
along outflows from very young stars. Kaufman & Neufeld (1996) predicted that between 
40 and 70% of the total shock luminosity is emitted in H 2 lines for shocks with pre-shock 
density lower than 10 5 cm -3 and shock velocities larger than 20 kms -1 , the other main 
contributions being in CO and H 2 rotational emission. This has been also observationally 
tested by Giannini et al. (2001) who measured the relative contribution of the different 
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species to the outflow cooling in a sample of class objects observed with ISO-LWS. 

We will discuss here the role of the H 2 cooling in the global radiated energy of the 
LI 157 outflow. From the best fit model obtained for the knots A and C, we have derived the 
total, extinction corrected, H 2 luminosity by integrating over all the ro- vibrational transitions 
considered by our model. L H2 is found to be 8.4xlCT 2 and 3.7xlCT 2 L for the A and C 
knots, respectively. Out of this total luminosity, the contribution of only the rotational lines 
is 5.6xlCT 2 (A) and 2.7xlO~ 2 (C) L , which means that in both cases they represent about 
70% of the total H 2 luminosity. 

N09 have found that the total luminosity of the H 2 rotational lines from S(0) to S(7), 
integrated over the entire L1157 outflow, amount to 0.15 L Q . If we take into account an 
additional 30% of contribution from the v>0 vibrational levels, we estimate a total H 2 
luminosity of 0.21 L Q . This is a 30% larger than the total H 2 luminosity estimated by 
Caratti o Garatti (2006) in this outflow, assuming a single component gas at temperature 
between 2000 and 3000 K that fit the NIR H 2 lines. 

If we separately compute the H 2 luminosity in the two outflow lobes, we derive Lh 2 = 
8.5xl0~ 2 L in the blue lobe and 1.3 xl0 _1 L Q in the red lobe. Comparing these numbers 
with those derived in the individual A and C knots, we note that the A knot alone contributes 
to most of the H 2 luminosity in the blue lobe. By contrast, the H 2 luminosity of the red lobe 
is distributed among several peaks of similar values. This might suggest that most of the 
energy carried out by the blue-shifted jet is released when the leading bow-shock encounters 
a density enhancement at the position of the A knot. On the other hand, the red-shifted 
gas flows more freely without large density discontinuities, and the corresponding shocks are 
internal bow-shocks, all with similar luminosities. 

The integrated luminosity radiated by CO, H 2 and OI in LI 157 has been estimated, 
through ISO and recent Herschel observations, as ~0.2 L Q (Giannini et al. 2001, Nisini et 
al. 2010), which means that H 2 alone contributes about 50% of the total luminosity radiated 
by the outflow. Including all contributions, the total shock cooling along the LI 157 outflow 
amounts to about 0.4 L , i.e. L coo ;/Lb D i ~ 5xl0~ 2 , assuming 1^=8.4 L for L1157-mm 
(Froebrich 2005). This ratio is consistent with the range of values derived from other class 
sources from ISO observations (Nisini et al. 2002). 

The total kinetic energy of the LI 157 molecular outflow estimated by Bachiller et al. 
(2001) amounts to 0.2 L without any correction for the outflow inclination angle, or to 
1.2 L if an inclination angle of 80 degrees is assumed. Considering that the derivation of 
the Lkm value has normally an uncertainty of a factor of five (Downes & Cabrit 2007), we 
conclude that the mechanical energy flux into the shock, estimated as L coo i, is comparable 
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to the kinetic energy of the swept-out outflow and thus that the shocks giving rise to the H 2 
emission have enough power to accelerate the molecular outflow. 

The total shock cooling derived above can be also used to infer the momentum flux 
through the shock, i.e. P = 2L coo i/V s , where V s is the shock velocity that we can assume, 
on the basis of the discussion in the previous section, to be of the order of 20 kms -1 . 
Computing the momentum flux separately for the blue and red outflow lobe, we derive P re d ~ 
1.7xl0~ 4 and Pu U e ~ l.lxlO~ 4 M Q yr" 1 kms -1 . In this calculation, we have assumed that 
the contribution from cooling species different from H 2 , as estimated by ISO and Herschel, 
is distributed among the two lobes in proportion to the H2 luminosity. If we assume that the 
molecular outflow is accelerated at the shock front through momentum conservation, then 
the above derived momentum flux should results comparable to the thrust of the outflow, 
derived from the mass, velocity and age measured through CO observations. The momentum 
flux measured in this way by Bachiller et al. (2001) is l.lxl0~ 4 and 2xl0 -4 M Q yr -1 kms _1 in 
the blue and red lobes, respectively, i.e. comparable to our derived values. It is interesting to 
note that the P determination from the shock luminosity confirms the asymmetry between 
the momentum fluxes derived in two lobes. As shown by Bachiller et al. (2001), the L1157 
red lobe has a 30% smaller mass with respect to the blue lobe, but a higher momentum 
flux due to the larger flow velocity. The northern red lobe is in fact more extended than the 
southern lobe: however, given the higher velocity of the red-shifted gas, the mean kinematical 
ages of the two lobes is very similar. 



5. Conclusions 

We have analysed the H 2 pure rotational line emission, from S(0) to S(7), along the 
outflow driven by the L1157-mm protostar, mapped with the Spitzer - IRS instrument. The 
data have been analysed assuming a gas temperature stratification where the H 2 column 
density varies as T _/3 and 2D maps of the H 2 column density, ortho-to-para ratio (OPR) and 
temperature spectral index (3 have been constructed. Further constraints on the physical 
conditions of the shocked gas have been derived in two bright emission knots by combining 
the Spitzer observations with near-IR data of H 2 ro-vibrational emission. Finally, the global 
H 2 radiated energy of the outflow has been discussed in comparison with the energy budget 
of the associated CO outflow. 

The main conclusions derived by our analysis are the following: 

• H 2 transitions with Ji OW er < 2 follows the morphology of the CO molecular outflow, with 
peaks correlated with individual CO clumps and more diffuse emission that delineates 
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the CO cavities created by the precessing jet. Lines with higher J are localized on 
the shocked peaks, presenting a morphology similar to that of the H 2 2.12/xm ro- 
vibrational emission. 

• Significant variations of the derived parameters are observed along the flow. The H 2 
column density ranges between 5xl0 19 and 3xl0 20 cm -2 : the highest values are found 
in the blue-shifted lobe, suggesting that here the outflowing gas is compressed due to 
the impact with a high density region. Gas components in a wide range of temperature 
values, from ~ 250 to ~ 1500 K contribute to the H 2 emission along individual lines 
of sight. The largest range of temperature variations is derived towards the intensity 
peaks closer to the driving source, while a more uniform temperature distribution, with 
T between 400 and 1000 K, is found at the tip of the northern outflow lobe. 

• The OPR is in general lower than the equilibrium value at high temperatures and spans 
a range from ~0.6 to 2.8, with the lowest values found in low temperature plateau 
regions between consecutive intensity peaks. As in previous studies, we also found the 
presence of regions at low OPR (1.5-1.8) but with relatively high temperatures. These 
might represent zones subject to recent shocks where the OPR has not had time yet 
to reach the equilibrium value. 

• Additional shock parameters have been derived in the two bright near-IR knots A 
and C, located in the blue- and red-shifted outflow lobes, where the mid- and near- 
IR H 2 data have been combined. The ratio between mid- and near-IR lines is very 
sensitive to the molecular plus atomic hydrogen particle density. A high abundance of 
atomic hydrogen (H/H 2 ~ 0.1-0.3) is implied by the the observed H 2 column densities 
if we assume n(H 2 ) values as derived by independent mm observations. With this 
assumption, the cooling lengths of the shock result of the order of 7xl0 14 and 10 15 cm 
for the A and C knot, respectively. The distribution of OPR values as a function of 
temperature and the derived abundance of atomic hydrogen, implies that the shock 
passing time is of the order of 100 yr for knot A and 1000 yr for knot C, given the 
assumption that the para-to-ortho conversion occurs through reactive collisions with 
atomic hydrogen. We find that planar shock models, either of C- or J-type, are not 
able to consistently reproduce all the physical parameters derived from our analysis of 
the H 2 emission. 

• Globally, H 2 emission contributes to about 50% of the total shock radiated energy in the 
LI 157 outflow. We find that the momentum flux through the shocks derived from the 
radiated luminosity is comparable to the thrust of the associated molecular outflow, 
supporting a scenario where the working surface of the shocks drives the molecular 
outflow. 



-17- 



This work is based on observations made with the Spitzer Space Telescope, which is op- 
erated by the Jet Propulsion Laboratory, California Institute of Technology under a contract 
with NASA. Financial support from contract ASI 1/016/07/0 is acknowledged. 

REFERENCES 

Arce, H. C, Santiago-Garcia, J., J0rgensen, J. K., Tafalla, M., & Bachiller, R. 2008, ApJ, 
681, L21 

Bachiller, R., & Perez Gutierrez, M. 1997, ApJ, 487, L93 

Bachiller, R., Perez Gutierrez, M., Kumar, M. S. N., & Tafalla, M. 2001, A&A, 372, 899 

Benedettini, M., Viti, S., Codella, C, Bachiller, R., Gueth, F., Beltran, M. T., Dutrey, A., 
& Guilloteau, S. 2007, MNRAS, 381, 1127 

Cabrit, S., et al. 1999, The Universe as Seen by ISO, 427, 449 

Caratti o Garatti, A., Giannini, T., Nisini, B., & Lorenzetti, D. 2006, A&A, 449, 1077 
Codella, C, et al. 2009, A&A, 507, L25 
Codella, C, et al. 2010, A&A, 518, L112 
Davis, C..J. & Eisloffel, J. 1995, A&A, 300, 851 

Dionatos, O., Nisini, B., Cabrit, S., Kristensen, L., Pineau Des Forets, G. 2010, A&A, 
submitted 

Downes, T. P., & Cabrit, S. 2007, A&A, 471, 873 

Eisloffel, J., Smith, M. D., & Davis, C. J. 2000, A&A, 359, 1147 

Flower, D. R., Le Bourlot, J., Pineau des Forets, G., & Cabrit, S. 2003, MNRAS, 341, 70 
Froebrich, D. 2005, ApJS, 156, 169 

Giannini, T., Nisini, B., & Lorenzetti, D. 2001, ApJ, 555, 40 

Giannini, T., McCoey, C, Caratti o Garatti, A., Nisini, B., Lorenzetti, D., & Flower, D. R. 
2004, A&A, 419, 999 

Giannini, T., McCoey, C, Nisini, B., Cabrit, S., Caratti o Garatti, A., Calzoletti, L., & 
Flower, D. R. 2006, A&A, 459, 821 



Gueth, F., Guilloteau, S., & Bachiller, R. 1996, A&A, 307, 891 
Gueth, F., Guilloteau, S., & Bachiller, R. 1998, A&A, 333, 287 

Gusdorf, A., Pineau Des Forets, G., Cabrit, S., & Flower, D. R. 2008, A&A, 490, 695 

Hirano, N., & Taniguchi, Y. 2001, ApJ, 550, L219 

Houck, J. R., et al. 2004, ApJS, 154, 18 

Kaufman, M. J., & Neufeld, D. A. 1996, ApJ, 456, 611 

Le Bourlot, J., Pineau des Forets, G., & Flower, D. R. 1999, MNRAS, 305, 802 

Lefloch, B., Cernicharo, J., Cabrit, S., Noriega-Crespo, A., Moro-Martin, A., & Cesarsky, 
D. 2003, ApJ, 590, L41 

Lefloch, B., et al. 2010, A&A, 518, L113 

Looney, L. W., Tobin, J. J., & Kwon, W. 2007, ApJ, 670, L131 
Maret, S., et al. 2009, ApJ, 698, 1244 
Molinari, S., et al. 2000, ApJ, 538, 698 

Neufeld, D. A., Melnick, G. J., & Harwit, M. 1998, ApJ, 506, L75 
Neufeld, D. A., Melnick, G. J., Sonnentrucker, P. et al. 2006, ApJ, 649, 816 
Neufeld, D. A., & Yuan, Y. 2008, ApJ, 678, 974 

Neufeld, D. A., Nisini B., Giannini T., et al. 2009, ApJ, 706, 170, N09 

Nisini, B., Benedettini, M., Giannini, T., Codella, C, Lorenzetti, D., di Giorgio, A. M., & 
Richer, J. S. 2000, A&A, 360, 297 

Nisini, B., Giannini, T., & Lorenzetti, D. 2002, ApJ, 574, 246 

Nisini, B., Codella, C, Giannini, T., Santiago Garcia, J., Richer, J. S., Bachiller, R., & 
Tafalla, M. 2007, A&A, 462, 163 

Nisini, B., et al. 2010, A&A, 518, L120 

Rodriguez-Fernandez, N. J., Tafalla, M., Gueth, F., & Bachiller, R. 2010, A&A, 516, A98 
Schofield, K. 1967, Planet. Space Sci., 15, 643 



- 19 - 

Wilgenbus, D., Cabrit, S., Pineau des Forets, G., & Flower, D. R. 2000, A&A, 356, 1010 

Rieke, G. H. & Lebofsky,, M. J., 1985, ApJ, 288, 618 

Zhang, Q., Ho, P. T. P., & Wright, M. C. H. 2000, AJ, 119, 1345 



This preprint was prepared with the A AS IATgX macros v5.2. 



-20 - 



Table 1. H 2 line brightness in selected positions 
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intensity averaged within a 20" FWHM Gaussian aperture around the considered position. Absolute 
errors are of the order of 20%. 

b Offsets in arcsec with respect to the LI 157- mm source. These positions are associated to individual 
near-IR/mm emission clumps, as indicated in Table 2. 

c 2cr upper limit 
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Table 2. Parameters derived from rotational diagrams' 
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a Typical uncertainty in the parameters is as follows: N(H 2 ) ±0.1 dex, OPR 
± 0.2, T± 50 K, £±0.2. 

h T co id and T warm are the temperatures derived from linearly fitting the S(0)- 
S(l)-S(2) and S(5)-S(6)-S(7) lines. T me d is the average temperature derived 
from a liner fit through all the H 2 lines. 

Association with mm clumps and with NIR-knots, as defined in Bachiller 
et al. (2001) and Davis & Eisloffel (1995), respectively. 
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Table 3. Inferred shock parameters' 1 



Parameter 
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a See text for details on the various parameters. 

A Model that assumes the H2 particle density 
from Nisini et al. (2007) and a variable H/H2 
ratio. 

B Model that assumes negligible hydrogen in 
atomic form. 
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Fig. 1. — L1157 maps of the 0-0 S(l)17.035 /xm and S(2)12.29 /zm lines observed with 
the Spitzer IRS-SH module. In both images contours are displayed in steps of 10% of the 
peak intensity value, which is 8.2xl0 -5 and 1.9xl0 -4 ergs _1 cm _2 sr _1 for the S(l) and S(2) 
maps, respectively. The average rms noise in the two images is of the order of 1.5 xlO -6 



ergs 'cm 2 sr 1 . 
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RA Offset (arcsec) 

Fig. 2. — (left) Map of the 1-0 S(5) line at 6.9 fim with superposed contours of the integrated 
CO (2-1) emission (from Bachiller et al. 2001). The position of the molecular condensations 
(bullets) identified in (Bachiller et al. 2001) are labelled, (right) Map of the 1-0 S(l) line at 
2.12 /im (Caratti o Garatti et al. 2006). The main group of NIR knots, labelled from A to D 
according to Davis & Eisloffel (1995) , are indicated. 
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Fig. 3. — Map of the total H2 column density. Values have been derived only in the pixels 
where at least four transitions with an S/N larger than 3 have been detected. Contours refer 
to the intensity of the S(l) line. 
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Fig. 4. — As in Fig. 3 for the map of the ortho-to-para ratio. 
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Fig. 5. — Maps of the fitted temperature power law spectral index (3 (a), and of the "cold" 
and "warm" temperature components, derived from a linear fit through the OPR corrected 
SfflVSM VSf2") lines and /Sffil /S>(7) lines, resnectivelv. Contours of the intensitv of the 
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Fig. 6. — H 2 rotational diagrams obtained at the offset positions indicated in each panel. The 
filled circles correspond to the observations while the solid line is the best LTE fit obtained 
in each position. The parameters of each fit are summarized in Tab. [2] 
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Fig. 7. — H2 rotational diagram showing the results of the NLTE model at different H2 
densities for the Spitzer observed rotational lines and three bright NIR lines at excitation 
temperatures between 6000 and 8000 K. The model assumes a slab where the temperature 
decreases as a power low with index /3 and where the o/p ratio varies as a function of 
temperature. The gas is assumed fully molecular (H/H 2 =0). Filled squares correspond 
to the observations in a 20" aperture centered towards the knot A in the south-east lobe 
(see FigJT]). Different symbols indicate the model results at different densities: the other 
parameters are indicated in the figure and are kept fixed. It should be noted how the near- 
IR lines are much more sensitive to density variations than the pure rotational lines, due to 
their different critical densities. 
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Fig. 8. — H 2 rotational diagrams showing the best model fits (open circles) through both 
the Spitzer and near-IR data (filled squares) in 20" apertures centered towards the knots A 
and C (see Fig ED). The gas is assumed fully molecular (H/H 2 =0). 
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Fig. 9. — As in Fig. |8] but keeping the H 2 density equal to the value of 4xl0 5 cm -3 as 
measured from SiO (Nisini et al. 2007), and varying the H/H 2 ratio. Agreement with the 
data is obtained with a neutral hydrogen abundance ~ 30% of the H 2 density. 



